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Abstract 


A lead-shielded HPGe detector and offline у-гау spectra of the residual product were used to 
measure the cross-section (CS) and ratios of isomeric CS (65/0) in "^Xe(n,2n)'?"*Xe reactions 
at different energies (13.5 MeV, 13.8 MeV, 14.1 MeV, 14.4 MeV, 14.8 MeV) relative to the 
?Nb(n,2n)""Nb reaction CS. The target was high-purity natural Xe gas under high pressure. 
The T(d,n)*He reaction produces neutrons. TALYS code (version 1.95) for nuclear reactions 
was used for calculations, with default parameters and nuclear level density models. The 
uncertainties in the measured CS data were thoroughly analyzed using the covariance analysis 
method. The results were compared with theoretical values, evaluation data, and previous 
experimental findings. CS data of the "^Xe(n,2n)'?"Xe and '**Xe(n,2n)'**8Xe reactions and the 
corresponding isomeric CS ratios at 13.5 MeV, 13.8 MeV, and 14.1 MeV neutron energies are 
reported for the first time. This research advances our knowledge of pre-equilibrium emission in 
the (n,2n) reaction channel by resolving inconsistencies in the Xe data. 
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(IRs); Covariance analysis 


1. Introduction 

Reliable knowledge of neutron-induced cross-sections (CSs) is required in nuclear power 
plants for energy and environmental protection, and in nuclear astrophysics reaction network 
calculations for vital input parameters [1-5]. Accurate measurements of CSs induced by neutrons 
at 14 MeV are desired for fusion reactors, nuclear transmutation, and nuclear model analyses 
[6-8]. Xe has nine stable isotopes [9]; 14-MeV neutrons can trigger nuclear reactions of Xe 
isotopes, including (n,a), (n,d), (n,p), (n,t), and (n,2n) [10,11]. The ^ Xe(n2n) "Хе CSs for 
such reactions were experimentally obtained by two laboratories [12,13] using sodium 
perxenate (Na4XeO¢-2H20) and quinol clathrate targets, respectively. In general, preparation of 
this compound is complicated; the composition is complex, and it is difficult to ensure purity 
(Ref. (12) reported an МаХеОв-2Н2О purity of 95.8%. Only one laboratory reported pure 
ground-state CSs in the "^Xe(n,2n)?*Xe reaction [12]. The excited-state, ground-state, total 
CS, and CS ratio for ^Xe(n,2n)' ?"Xe reactions correspond to a single energy value (14.4 
MeV or 14.6 MeV). Extensive research on the excitation functions for !2^126128.130.132.36y e(n, 
2n) and P9? ?Xe(n. p) reactions in the range of 13-15 MeV was recently reported by Bhatia 
et al. [14], Bhike et al. [15,16], and Luo et al. [17,18]. 

Our goal was to analyze four aspects of the '**Xe(n,2n) reaction CS. First, the target “Хе 
isotope belongs to a transitional region (N = 82, below the closed neutron shell); an even "*Xe 
nucleus for such shell closure is a good indicator of single-particle excitations [17]. Second, the 
half-life of the metastable state of the '**"Xe pair is shorter than that of the ground state. This 
pair exhibits isomeric transition (IT) decay (Fig. 1). Third, published data on the '**Xe(n,2n) 
reaction CS are scarce owing to difficulties in preparing gaseous targets for irradiation. Fourth, 
significant discrepancies have been observed in published results [12,13], mainly due to the 
following three factors: 

1. Systematic inconsistencies are caused by differences in experimental techniques and 
conditions (neutron field characteristics, target sample, neutron monitoring method). 
Previous studies have used sodium perxenate (Na4Xe0O6:2H20) and the solid quinol 
clathrate of Xe as the target materials [12, 13]. For neutron flux monitoring, Sigg and 
Kuroda (1976) [12] used a low-threshold reaction: "Al(n,à) Na (Е,=3.25 MeV). 


Kondaiah et al. (1968) [13] used "Xe(n2n) ^ *9Xe and P*Xe(n2n)?" "9 Xe reactions 
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iii. 


to monitor "^Xe(n,2n)?"Xe and "^Xe(n,2n) ? 9 Xe nuclear reactions (Table 1). 
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Figure 1. Level schemes of isomeric pairs of "Хе [9]. АП energies аге in keV. 


There is a significant shortage of decay data. The '“Хе(п,2п) "Хе reaction CSs [12, 13] 
exhibited characteristic y-ray intensities (14% and 13.55% at 233 keV, respectively 
(Table 1)). Recently, 10.12% was also reported [9]. The "^Xe(n,2n) "Xe reaction CSs 
were measured by Sigg and Kuroda (1976) [12], and Kondaiah et al. (1968) [13] using 
characteristic y-ray intensities of 35% and 35.7%, respectively, for the same 
characteristic y-ray of 81 keV (Table 1). However, 36.9% was obtained more recently 
[9]. As 79.61 keV and 80.998 keV y-ray energies are so close that they cannot be 
distinguished by the detector, the intensity of the 81-keV gamma ray should be higher 
(37.34% = 0.44% + 36.9%). The CS of '?"Xe radionuclide product with a 2.26-d half- 
life was also determined [12, 13], recently reported as 2.198 d [9]. 

There are differences in the total CS measurement methods. Only two methods currently 
exist for calculating the total CS. 1) The excited-state and pure ground-state CSs can be 
added after measuring them separately [12]. 2) The ground state can be measured after 
full decay of the excited state via IT, and the result can be considered as the total CS [13]. 
The latter method satisfies the conditions that the excited-state half-life is less than that 
of the ground state, preferably much less, and IT = 100%. If these conditions are not 


satisfied, it is not advisable to replace the total CS with the obtained ground-state CS. 


Analysis of the existing data led us to conclude that the isomeric ratios (IRs) and CSs for 
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the (n,2n) reaction of '**Xe in the 13-15 MeV range must be further analyzed. We used residual 
nuclear decay analysis [19—21] and natural Xe with high-purity pressure as the target to measure 
the "^Xe(n2n)?"Xe, U*Xe(n2n)?*Xe, and "Хе(п,2п) Хе CSs reactions and their 
corresponding isomeric-state CSs in the 13-15 MeV range. The results were compared with 
theoretical values obtained using TALYS-1.95 [22], evaluation data, and previous experimental 


results. 


Table 1. Literature indicating experimentally obtained CSs for ?^Xe(n,2n) "Хе reactions 


En 
Reaction Decay data Sample Detector Monitor reaction (MeV) с (mb) Ref. 
e 
^X e(n2n) 37. T1n=2.26 d, Ey-232.8 К Na4XeO«2H GeLi ?7A1(n,a^^N 12 
e(n,2n) e р 6 d, Еу=232.8 keV, a4XeOs2H;0 eLi (n,a) Na 14.6 655460 906483" [12] 
Iy-14 96 
7,572.26 d, Ey-233 keV, inol clathrate Geli "еХе(п,2п) ^" 9x 13 
р ‚ Еу eV, Quinol clathrate eLi e(n,2n) e 144 665480 39041078 [13] 
113.55 % 
P^Xe(n2n)"*Xe 71=5.29 d, Ey-81.0 keV, y= NaiXeO&2H;O  GeLi 7 Al(n,a) Na 14.6 805490 755-84 [12] 
35% 14.4 1695+226° 1366+202° [13] 
P^Xe(nj2n)P?""?Xe 7,5952 Ey-81.0 keV, [= NasXeOo-2H GeLi 27A l(n,a.))* 12 
e(n,2n) e 575.29 d, Ey-81.0 keV, № a4XeOs2H;0 eLi (n,a) Na 14.6 1460-110 — 16615118 [12] 
35% 
Tin=5.2 Ey=81 k h= inol clathrate GeLi У e(n2n) "9x 13 
da 3 7 d, Ey=8 eV, [= Quinol clathrate GeLi e(n,2n) e i4 23604240 22564229" [13] 
* 0 


a: corrected results using characteristic y-ray intensity 
b: Reference 13 does not provide ground-state CS values. This value was obtained by subtracting the excited-state CS value from the total CS value. 


2. Experimental Methods 
2.1. Sample preparation 

Naturally pure Xe (99.999954, 2.3—5.2 g) pressurized to 120—250 atm was placed in а 1-mm 
thick spherical stainless steel container 20 mm in diameter. The samples and monitors (niobium 
foil: Ф = 20 mm, 0.12-mm thick, 99.99% pure; aluminum foil: Ф = 20 mm, 0.3-mm 
thick,99.999% pure), and zirconium foil (Ф = 20 mm, 0.2-mm thick, 99.95% pure) were 
combined as ZrNbAI-Xe-AINbZr for irradiation. The samples were mounted at different angles 
(0°, 45°, 90°, 110°, 135?) to the incident beam. The neutron-producing target edge was 5 cm 
from the sample center. Due to its high threshold energy, the ?Nb(n,2n)??"Nb (threshold energy 
8.972 MeV) reaction was used to monitor the corresponding threshold energies of 
*Xe(n,2n)??"£Xe reactions (Еһ = 8.853 MeV and 8.618 MeV). Low-threshold reactions such 
as (n,a) and (np) were monitored via a ”Al(n,a)” Na reaction (threshold energy = 3.249 MeV). 
The CS of the ?Nb(n,2n)"?"Nb reaction was obtained from the IRDFF-II data library [23]. To 
minimize the influence of the "?Xe(n,y)' ?"*Xe reaction on the "^Xe(n,2n)'?"Xe reaction, 
the samples were covered with Cd foil (1-mm thick, 99.9595 pure). The details of this approach 
are presented in our previous studies [17,18]. 

The neutron generator facility was provided by the China Academy of Engineering Physics 
(CAEP) for sample irradiation at (4-5)x10!° n/s yield (for 2 В). The d+*>H->n+*He+17.6 MeV 
reaction provided neutrons at 14 MeV with deuteron beam current and energy of 200 uA and 
134 keV, respectively. The tritium-titanium target thickness was 2.65 mg/cm?. A beam of a- 


particles at 135? was used to correct the neutron flux. 


2.2. Neutron energy and detection efficiency 


In the experiment, Zr and Nb were used to measure the CS ratios of ?Nb(n,2n)?"Nb and 
?7r(n,2n)?"" Zr, respectively, to determine the neutron energy. At irradiation angles of 135°, 
110°, 90°, 45°, and 0°, the energies were 13.5 MeV, 13.8 MeV, 14.1 MeV, 14.4 MeV, and 14.8 
MeV, respectively. For verification purposes, the neutron energies were determined using 
different methods [24]. Considering the sample distance from the target, size, and beam radius, 
the uncertainty in the neutron energy was 0.2 MeV [25]. The results obtained using the CS ratio 


method were consistent in terms of the experimental errors. The detector efficiency was 


calibrated using ‘Eu, "Ba, "Cs, and 2а standard sources. The Monte Carlo method was 
used to adjust the data for the geometrical differences between the calibration sources applied 
to the HPGe detector and the Xe sphere efficiency determination. The details are provided in 


our previous studies [17,18]. 


2.3. Radioactivity measurement 


Five samples (No.1—5) were allowed to decay after irradiation for 2.9 d, 3.9 d, 5.0 d, 4.0 d, 
and 5.0 d, respectively. The y-ray activities of "Хе, !?*Xe, and ?"Nb nuclei were measured 
with a coaxial HPGe detector, with crystal diameter, length, relative efficiencies, and energy 
resolution (1.332 MeV) of 70.1 mm, 72.3 mm, 6895, and 1.69 keV, respectively. Typical y-ray 
spectra obtained after 2.9 d of irradiation are shown in Fig. 2. A y-spectrum analysis was 


performed using ORTEC® GammaVision® software to estimate the peak area [26]. 
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Figure 2. Xe y-ray spectrum after 2.9 d of decay after irradiation. Data were collected for 23.4 
h. 


Table 2 displays the natural isotope abundance and product decay characteristics. 


Table 2. Details of current "^Xe(n,2n) nuclear reaction CS measurements (taken from ENSDF 
(2022) [9]). 


Isotope Reaction E- Mode Product Ру (кеу)  ly(96) 
Abundance (%) threshold decay (%) ^ half-life 
(MeV) 
10.43572; P^Xe(n2n)?"Xe 8.853 IT(100) 2.198 аз 233.22 10.125 
10.43572; P^Xe(n2n)?*Xe 8.618 p(100) 5.2475 4, 80.998 36.9; 
79.61 0.44 5 
100 9?Nb(n,2n)?"Nb 8.972 EC (100) 10.15 d? 934.44 99.15; 


Uncertainties are indicated by the lower index and italicized numbers (10.43572;% indicates 
(10.4357 + 0.0021)%, 2.198 аз indicates (2.198 = 0.013) d. 


3. CS calculations and associated uncertainties 
3.1. CS calculation 


The reaction CSs can be expressed as [17, 18] 


_[Sel,nKMD]y, [AAFC], 


om O Nb (1) 
[SeI,nKMD], [AAFC],, 


where Nb and x are the monitored and measured values, respectively; e is the full-energy peak 


efficiency; „Г, represents the gamma-ray intensity; 77 is the abundance of the target nuclide; M is the 


t —A(t *t;) T 


mass of the sample; D- e ^ —e represents the counting collection factor; S =1-е ^" is 
the growth factor of the product nuclide; Т denotes the irradiation time; t; is the cooling time; #2 is 
the measurement time; А is the atomic weight; C represents the measured full-energy peak area; 4 
is the decay constant of the residual nucleus; K is the correction factor for decaying nuclei during 


irradiation time (T), which is divided into small time intervals Az;. This factor is calculated as 


K= [»: & — e ^^ )е ^^ ]/ (SQ) , where iis the mean neutron fluence rate of the i part, and is 
m 


regarded as a constant (when и is sufficiently large). T; is the time from the end of the i part to the 


n 
end of total irradiation; ф= At, / Т is the mean neutron fluence rate within the irradiation 
i=l 

time. F is the total correction factor, expressed as 


Е=ЕхЕХЕ, (2), 


where Ро, Fs, and F, are the correction factors for the solid angle of the neutron flux, self- 
absorption of the specific y-ray energy, апа sample-counting geometry, respectively. Fs, Fg, and 


Ко were estimated using the characteristics of spherical samples, as described in Section 3.1 of 
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Ref. [18]. Table 3 presents the self-absorption correction factors. 


Table 3. Correction variables for sample self-absorption at different y-ray energies 


Gamma-ray ulp (cm?/g) u(cm^) Sample — 
pa Fe Xe Fe | No Хе por factor fs 
80.998 0.5840 3.5521 4.5906 1 4.4117 1.2420 8.288 

2 2.0804 0.5857 4.103 
3 4.2326 1.1916 7.956 
4 1.9529 0.5498 3.889 
3 1.9369 0.5453 3.863 
233.22 0.1340 0.3108 1.0533 1 0.3860 1.2420 1.390 
2 0.1820 0.5857 1.204 
3 0.3703 1.1916 1.375 
4 0.1709 0.5498 1.195 
5 0.1695 0.5453 1.194 


To obtain the pure ground-state CSs for "*Xe(n,2n)'?*Xe reactions, С, in Equation (1) 
was set to а peak area of 80.998 keV minus the contribution from Хе via 


133m e 00091338 Хе (counting Ci). Ci can be expressed using the decay of an artificial 


nuclide, expressed as [18] 


2 2 
C = Eee eu (A, SEEN -ASe D) (3) 
0 (&7AJS D LEAK, 


mmm g 


where subscripts g and m are the ground and metastable states, respectively; Ри» is the 


percentage of metastable-state disintegrations that produce ground-state nuclides; C, is the 


-å 


energy peak area of the metastable state; 5, —1—e ^" and $, =1- e ^ ; Im and I, are the 


intensities of y-rays, where =„ and & are the full-energy peak efficiencies of the metastable and 
ground states, respectively; K, is the neutron-injection fluctuation indicator of the metastable 


state; Fm and Е. denote the correction factors of the metastable and ground states, respectively; 


4,4 —An (4 +t) 27% (fj +t) 


Dm and Dg are coefficients defined as D, —e ^"^ —e and D, =e ^^ - 


3.2. Associated experimental uncertainties 


When calculating the main uncertainties of CS determination, we assumed that only the 
half-life uncertainty contributed to the uncertainty in the timing factor. The following 
parameters were used for this purpose: C. as the y-ray counting statistics, М, уь as the target 
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masses, /. м as the relative y-ray intensities, 74, as the target isotopic abundance, ow, as the 
standard CS, & ль as the efficiencies, and 5, уь, and D, as timing factors. The subscripts Nb 
and x correspond to the monitored and measured reaction-related terms, respectively. 

For generative nuclei "Хе and Хе, the uncertainty of correction factor F is estimated 
as 2.096 and 3.095, respectively. These uncertainties were applied to assess the total uncertainty 
of the experimental CSs [27-30]. Table 4 and Table 5 highlight the uncertainties in the 
parameters contributing to the reaction CS values, which were used to extend the covariance 
matrix between different energy levels. All samples were investigated using the same setup and 
standard CS reaction. Thus, the detector efficiency and standard CS accuracy were consistent 
for all neutron energy correlations. After calculation of fractional uncertainties, covariance 
analysis calculated the correlation coefficients between each energy-related property. The 
coefficients and corresponding energies are presented in Table 6. The CS covariance matrix 
and total uncertainty were generated for energy-level pairs (e.g., (Oxi,0,j)) using the data from 


Tables 4—6 by adding 15 subset matrices in Equation 4 [29]: 


Соу(с 


хі? 


a) =}, 2, Ax, x Cor(Ax,, Ax) x Ах, (4). 
i j 


The generated covariance matrix was [5x5]. The uncertainty in the experimental CS value 


was determined as (Equation 5) [29] 
Cov(o.,,0.,) = (A0)? (5). 


From the covariance matrix and total uncertainty between neutron energies, the correlation 
matrix [5x5] was derived using Equation 6 [29]: 


Cov(O iy) (6). 


Cor(o 
(Aa (Ao, 


o0) = 


The experimentally obtained reaction CS values and their correlation matrices and uncertainties 
are presented in Table 7 and Table 8. 
Table 9 shows experimental values for !3Хе(п,2п)! "Хе reaction CSs and their IRs, and the 


monitor reaction CSs at different neutron energies. 
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Table 4. Uncertainties (%) in parameters for ^Xe(n,2n)?"Xe reaction CSs at 13.5-14.8 MeV 


Fractional uncertainties (96) 


Attributes spo ЕЕЕ 

(x) 13.5 MeV 13.8 MeV 14.1 MeV 14.4 MeV 14.8 MeV 
(Axi) (Axj) (Axx) (Axi) (Axn) 
Cy 0.1236 0.1471 0.1349 0.1750 0.0726 
Cub 0.8523 0.8102 0.8083 0.8164 0.9315 
Ix 1.4822 1.4822 1.4822 1.4822 1.4822 
Inb 0.0403 0.0403 0.0403 0.0403 0.0403 
Mx 0.0401 0.0869 0.0816 0.0876 0.0385 
Мм 0.0111 0.0109 0.0112 0.0110 0.0112 
Nx 0.0201 0.0201 0.0201 0.0201 0.0201 
ONb 0.6327 0.5615 0.5461 0.5427 0.5645 
& 2.7201 2.7201 2.7201 2.7201 2.7201 
ЕМЬ 2.5038 2.5038 2.5038 2.5038 2.5038 
Sx 0.5836 0.5836 0.5836 0.5836 0.5836 
Sub 0.1965 0.1965 0.1965 0.1965 0.1965 
Dx 0.4188 0.4197 0.2302 0.2335 0.0419 
Dw 0.1857 0.1857 0.1858 0.1858 0.1860 
Fy 2.0000 2.0000 2.0000 2.0000 2.0000 

Total error (%) 4.65 4.63 4.62 4.62 4.63 


Table 5. Uncertainties (%) in parameters for ?*Xe(n,2n) Хе reaction CSs at 13.5-14.8 MeV 


. Fractional uncertainties (94 
Attributes i) 


(x) 13.5 MeV 13.8 MeV 14.1 Мем 14.4 MeV 14.8 Меу 
(Axi) (Ах) (Axi) (Axi) (Axa) 
Cy 0.1198 0.1224 0.1115 0.1174 0.0936 
Смь 0.8523 0.8102 0.8083 0.8164 0.9315 
1, 1.2854 1.2854 1.2854 1.2854 1.2854 
Тм 0.0403 0.0403 0.0403 0.0403 0.0403 
Mx 0.0201 0.0869 0.0816 0.0876 0.0385 
Мм 0.0111 0.0109 0.0112 0.0110 0.0112 
Nx 0.0201 0.0201 0.0201 0.0201 0.0201 
ONb 0.6327 0.5615 0.5461 0.5427 0.5645 
& 3.9743 3.9743 3.9743 3.9743 3.9743 
£Nb 2.5038 2.5038 2.5038 2.5038 2.5038 
Sx 0.0094 0.0094 0.0094 0.0094 0.0094 
SNb 0.1965 0.1965 0.1965 0.1965 0.1965 
Dx 0.0027 0.0039 0.0039 0.0039 0.0052 
Dwp 0.1857 0.1857 0.1858 0.1858 0.1860 
Ех 3.0000 3.0000 3.0000 3.0000 3.0000 


Total error (%) 5.83 5.82 5.82 5.82 5.83 


Table 6. Correlation coefficients of parameters at 13.5-14.8 MeV for "^Xe(n2n) "Хе 
reaction CS. The Zy, Inb, Их, Омь, £v», Sno, and Оль subsets are fully correlated at these energies. 
Correlation coefficient (Ax, Ax) 


Cor(Ax, Ax) 

Cy Cp Mx Мм & Sx Dx Р, 
(Ax; Ах) l l 1 l 1 1 1 1 
(Axi, Axj) 0 0 0 0 0.9987 0 0 0 
(Ax;, Axx) 0 0 0 0 0.9987 0 0 0 
(Axi, Ах) 0 0 0 0 0.9987 0 0 0 
(Axi, Ах») 0 0 0 0 0.9987 0 0 0 
(Ах, Ах) | | | | | 1 1 | 
(Ах, Axx) 0 0 0 0 0.9987 0 0 0 
(Ax;, Axi) 0 0 0 0 0.9987 0 0 0 
(Ax;, Axn) 0 0 0 0 0.9987 0 0 0 
(Axx, Ax) ] ] ] 1 1 1 1 1 
(Axx, Axi) 0 0 0 0 0.9987 0 0 0 
(Ахь, Axn) 0 0 0 0 0.9987 0 0 0 
(Ax), Ах) 1 | | 1 | | | | 
(Axı, Ах») 0 0 0 0 0.9987 0 0 0 
(Axn, Ах,) 1 1 1 | 1 | | | 


Table 7. Experimental CSs (mb) of '**Xe(n,2n)'**"Xe reaction and their overall uncertainties 
and correlation matrices 


Neutron energy Cross-section Aox г . 
Correlation matrix 
En (MeV) сх (mb) (90) 
13.5=0.2 884-41 4.65 1.0000 
13.8=0.2 888-41 4.63 0.7565 1.0000 
14.1+0.2 909+42 4.62 0.7586 0.7593 1.0000 
14.4+0.2 916+42 4.62 0.7581 0.7588 0.7610 1.0000 
14.8+0.2 928+43 4.63 0.7562 0.7569 0.7591 0.7586 1.0000 


Table 8. Experimental CSs (mb) of "^Xe(n,2n)'?*Xe reaction and their overall uncertainties 
and correlation matrices 


Neutron energy Cross-section Aox . . 
Correlation matrix 


En (MeV) Ox (mb) (76) 

13.5=0.2 824-48 5.83 1.0000 

13.8=0.2 810447 5.82 0.7125 1.0000 

14.1+0.2 786-46 5.82 0.7125 0.7129 1.0000 

14.4+0.2 767+45 5.82 0.7123 0.7127 0.7127 1.0000 
14.8+0.2 747444 5.83 0.7105 0.7108 0.7108 0.7106 1.0000 


Table 9. Experimental CSs (mb) of ?^Xe(n,2n) ?""*Xe reaction and their IRs 


| Monitor reaction 
This study às 3 
En (MeV) Nb(n,2n) "Nb 
CS (mb) IRs (65/04) CS (mb) Reference 

13.5+0.2 1708-63 1.07+0.08 453.0+2.9 [24] 
13.80.2 1698-63 1.10+0.08 457.4+2.6 [24] 
14.1+0.2 1696-62 1.16+0.09 459.7+2.5 [24] 
14.4+0.2 1683-62 1.20+0.09 460.1+2.5 [24] 
14.8=0.2 167561 1.24+0.09 460.2+2.6 [24] 


4. Theoretical calculations 


TALYS software is widely used to analyze pre-equilibrium mechanisms, direct nuclear 
reactions, and compound nuclei (CN). It can evaluate nuclear reactions in the range of 0.001— 
200 MeV for target nuclei with masses of 12 and higher. TALYS also calculates the nuclear 
reaction CSs. Theoretical calculations for different ratios of ?^Xe(n,2n) "Хе reactions 
(O»/og) were performed for different ranges of neutron energies using TALYS-1.95 (ver. 1.95) 
[22]. We used the default parameters for our calculations, except for six distinct nuclear level 
density models (NLDs), which were divided into two groups: three microscopic levels and three 
phenomenological densities. The phenomenological level densities included: 1) Idmodel 1, with 
a constant temperature and Fermi gas model, and excitation energy divided into higher and 
lower energy regions in the Fermi gas model and with constant temperature laws, respectively; 
2) Idmodel 2 is based on the back-shifted Fermi gas model, which uses the Fermi gas expression 
across the region; 3) Idmodel 3 is a generalized superfluid model that considers superconductive 
correlations based on Barden-Cooper-Schrieffer theory [31-32]. The microscopic level 
densities are 1) ldmodel 4 is made up of Goriely's tables derived from the Skyrme force 
(microscopic level densities) [33]; 11) Idmodel 5 consists of Hilaire's combinatorial tables 
derived from the Skyrme force [33]; iii) ldmodel 6 is the Gogny force, derived from Hilaire’s 


combinatorial tables [33]. Examining the spin distribution is feasible using relative feeding as 
a ‘probe’. o2(E,) indicates the angular momentum distribution of the level density, 


expressed as 


A? 
67 (E,) = 0.01389 —./a(E, — ^) (7), 
а 
where A indicates the mass number; Ё; is the true excitation energy, A represents the empirical 
parameter and energy shift, which are nearly equal to the pairing energy of the observed nucleus 
odd-even effects; о is the energy-dependent level density parameter accounting for damping at 
higher excitation energies and shell effects at lower energies. Without shell effects, the 


asymptotic level density parameter ( à ) is equal to æ [34]. 


5. Results and discussion 


Tables 7-9 show the CSs of the "*Xe(n,2n) ?"Xe, ?^Xe(n2n) "Хе, and "^Xe(n2n) Хе 
reactions and the corresponding isomeric-state CS ratios obtained using the activation method. 
The data of the excited state and pure ground state were compared to data from the literature 
[12,13] and data obtained by TALYS-1.95 at six energy-level densities [22]. The 
P*Xe(n2n) Хе CS was also compared to data from the ЕМОЕ/В-УШ.О [35], JEFF-3.3 [36], 


JENDL-4.0 [37], ROSFOND [38], and CENDL-3 [39] databases. 


5.1. ?^Xe(n,2n)?"Xe reaction 


Two studies [12,13] have reported "^Xe(n,2n) experiments for the "Хе reaction. 
However, they provide the CS at only one energy value (14.4 MeV or 14.6 MeV). For a y-ray 
with 233.22 keV energy, Sigg and Kuroda [12] used 7,7 14%, and Kondaiah et al. [13] used 
I,— 13.6%. However, other studies indicate that data measurement accuracy was higher with 
[,= 10.12%. In this study, y-ray energy of 233.22 keV with an intensity of /,= 10.12% emitted 
during the decay of "Хе was used to obtain the '“*Xe(n,2n)'**"Xe reaction CS. Thus, the CS 


values were modified according to the following equation: 


I 
Ве 0 (8), 


corrected ^ — I literature 
y (new) 


where Ojiteraure is the CS from the literature; Г (ола) is the y-ray intensity at 233.22 keV, and Lynew) 
is 10.12%. The experimental data are strongly correlated with the corrected experimental results 
[12, 13], as shown in Fig. 3. However, the data are below the TALYS-1.95 excitation curves 
(Idmodels 1-6) for neutrons with energies of 13.8—14.8 MeV. These findings were first 


published at 13.5-MeV, 13.8-MeV, and 14.1-MeV neutron energies. 
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Figure 3. ?^Xe(n,2n) "Хе reaction CSs 

5.2. ?^Xe(n,2n)?**Xe reaction 

A thorough literature review yielded only one measurement of the "^Xe(n,2n) ?*Xe 
reaction using y-ray energy of 81 keV (4= 35%) [12]. However, many investigations have used 
1,= 36.9% for gamma-rays with the same energy. Production of ће Хе nucleus resulted in 
two characteristic gamma-rays with similar energies, 79.61 keV (1, — 0.44%) and 80.998 keV 
(4,7 36.9%). Our detectors could not resolve these energies; the intensity should be equal to 
their sum (37.3496). Accordingly, Equation (8) was used to modify the CS values from Ref. 
[12]. During the "Хе decay, the emitted 80.998-keV y-ray was used to calculate the 
P^Xe(n,2n) Хе reaction CSs after subtracting the contribution of the "^Xe(n,2n)'?"Xe 
reaction via the isomeric transition (IT). Almost 50% of the 80.998-keV у-гау counts were from 
the isomeric transitions. Figure 4 shows the theoretical assessments of the excitation functions 
and the experimental results (Ref. [13] does not provide ground-state CS values, which are 
equal to the total CS value minus the excited-state CS value.). The measured values exceeded 
those calculated using the TALYS-1.95 software with ldmodels 1-6. However, our 
experimental and theoretical results show the same dependence on neutron energy. At 14.6 
MeV, our measurements agreed with the corrected values obtained by Sigg and Kuroda [12]. 


At 14.4 MeV, the results in Ref. [13] significantly exceeded our experimental and simulated 


results (using TALYS-1.95, Idmodels 1-6). For this reaction, the CS values span a wide energy 


range. 
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Figure 4. "^Xe(n,2n) ?*Xe reaction CSs 

5.3. *Xe(n,2n)'?Xe reaction 

There are two published reports on the "*Xe(n,2n) Хе reaction at a neutron energy of 14 
MeV that can be used to validate theoretical calculations and experimental results. The shapes 
of the excitation curves obtained from the ENDF/B-VIII.0 [35], (JEFF-3.3 [36], CENDL-3 
[39]), JENDL-4.0 [37], and ROSFOND [38] databases were almost identical to those obtained 
by TALYS-1.95 [22] at ^15 MeV (Fig. 5). We conclude that the experimental data agree with 
the ROSFOND [38] results, TALYS-1.95 with Idmodel 1, and the corrected experimental 
results of Sigg and Kuroda [12]. At 14.4 MeV, the values [13] were 570 mb higher than ours. 
The evaluation excitation curve taken from JENDL-4.0 [37] is above the corrected results of 
Sigg and Kuroda [12], our experimental results, and results from the other databases 


[35,36,38,39]. Data for the 13.5—14.1 MeV neutron energy range are reported for the first time. 
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Figure 5. '**Xe(n,2n)'**Xe reaction CSs 

5.4. Isomeric CS ratio 

The ratios of the isomeric CSs were calculated using the measured "^Xe(n,2n)' ?"£Xe CSs. 
The CS ratios (с/с) for the isomeric "Xe pair created from the "*Xe(n,2n) reaction 
exposed to incident neutrons with energies of 13.5 + 0.2 MeV, 13.8 + 0.2 MeV, 14.1 + 0.2 MeV, 
14.4 € 0.2 MeV, and 14.8 + 0.2 MeV are 1.07 + 0.08, 1.10 + 0.08, 1.16 + 0.09, 1.20 + 0.09, and 
1.24 + 0.09, respectively (Table 9). Figure 6 shows the corrected isomeric CS ratios from the 
literature [12,13] and the excitation curves obtained using TALYS-1.95 calculations for 
different density-level models. The isomeric CS ratios increased with energy. In the 14-MeV 
region, the isomeric CS ratio is directly correlated with the neutron energy, indicating that 
higher excitation energies result in high-spin isomers (11/2.—3/2^) [40-43]. The obtained 65/6; 
values were consistent with the modified values from Sigg and Kuroda [12], but slightly lower 
than the values obtained using Ildmodels 1—6 and higher than the results of Ref. [13]. Our data 
for neutrons with 13.5—14.1 MeV energies are the first available in the literature to the best of 


our knowledge. 
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Figure 6. IRs of ?^Xe(n,2n) "Хе reactions as a function of neutron energy 

6. Summary and Conclusions 

The сь, Og, and Om+g for the PX e(n,2n) "Хе, ЗХ е(п,2п) ?*Xe, and P*Xe(n2n) ?Xe 
reactions and the corresponding isomeric CS ratios (Om/Og) were determined from the neutron 
activation and offline y-ray spectra of the residual products at 13.5-14.8 MeV using the 
?Nb(n2n)""Nb reaction as a reference. We thoroughly analyzed the uncertainties of the 
experimentally obtained CSs using covariance analysis. Experimental data were compared to 
data from the literature and CS values obtained using TALYS-1.95 software and from the 
ENDF/B-VIILO, JEFF-3.3, JENDL-4.00 ROSFOND, and CENDL-3 databases. The 
theoretically obtained CS values for the "^Xe(n,2n) reaction at six different NLDs (to generate 
nuclear excited and ground states) disagreed with the experimental results. However, the total 
CSs values were in agreement. The magnitude of the calculated CSs differed due to variations 
in the nuclear inputs and reaction models used for the simulations. The reaction mechanism 
involves photon and particle emissions in the CN. The CN CS is determined by the OPs, NLDs, 
and ySFs. Significant discrepancies were found in the experimental values; however, these were 
observed only for the pure ground state and total CSs as they are the most affected by the excited 
states. The CSs at neutron energies of 13.5—14.1 MeV are presented for the first time for the 


*Xe(n2n) "Xe, ^Xe(n 2n) Xe, and P^Xe(n2n) Хе reactions, significantly extending 
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the currently published energy ranges. These values can benefit nuclear databases for exploring 


fundamental nuclear applications and nuclear reaction models. 
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